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ABSTRACT 

Using the Long Wavelength Spectrometer on Keck, we have imaged the gravitationally Icnsed radio 
quiet quasi-stellar object (QSO) 2237-1-0305 at 8.9 and 11.7 /itm for the first time. The mid-infrared 
flux ratios are inconsistent with the optical flux ratios, but agree with the radio flux ratios and with 
some published gravitational lens models. These flux ratios indicate that the infrared emission is not 
affected by microlensing, which rules out the synchrotron emission model. The infrared emission is 
likely produced by hot dust extended on a length scale of more than 0.03 parsecs. The spectral energy 
distribution further implies a narrow range of dust temperatures, suggesting that the dust may be located 
in a shell extending between ^ 1 and 3 pc from the nucleus, and intercepting about half of the QSO 
luminosity. 

Subject headings: gravitational lensing — infrared: galaxies — quasars: QSO 2237+0305 



1. INTRODUCTION 

QSO 2237+0305 (z=1.695, Huchra et al. 1985), referred 
to hereafter as "the Einstein Cross," is the first gravi- 
tationally lensed QSO with confirmed variations due to 
lensing by the stars in the lens galaxy (z=0.0395), re- 
ferred to as "microlensing" (Irwin et al. 1989). Op- 
tical variations of ~ 0.5 magnitudes due to microlens- 
ing are observed at least once per year (Racine 1992), 
which implies the optical emission region is smaller than 
an Einstein radius, i?E = Fl x 10-^^M^''^/i~^/^cm, where 
M is the mass of the lensing star in solar masses and 
h = iJo/(100kms~^ Mpc~^). The peak magnification and 
time scale of a high amplification event (caused by a caus- 
tic crossing the QSO) is strongly dependent on the size of 
the QSO. An event in 1988 has been used to limit the opti- 
cal source size to about 10^'^ cm (Wambsganss, Schneider, 
& Paczyhski 1990, Ranch & Blandford 1991, Wyithe et al. 
2000). 

Typically QSOs have a change in spectral slope near 
1/im in the rest frame, which may imply that different 
physical components contribute longward and shortward 
of this wavelength (Neugebauer et al. 1987, Sanders et 
al. 1989, Barvainis 1993, Begelman 1994). One way to 
distinguish these components is reverberation mapping: 
observations of a few Seyfert galaxies indicate that the 
infrared emission region is more extended than the op- 
tical (Barvainis 1992, Nelson 1996). However, there ex- 
ist no such direct constraints on the size of the infrared 
emission in QSOs; recent observations indicate this re- 
gion may be rather compact, although the data are am- 
biguous (Neugebauer & Matthews 1999). Another tech- 
nique for constraining the size of the emission region in 
QSOs is gravitational microlensing. Sources larger than 
an Einstein radius have smaller fluctuations than sources 
smaller than an Einstein radius (Kayser, Refsdal, & Sta- 



bell 1986). Despite its potential for probing the QSO cen- 
tral engine, microlensing has not yet been looked for in 
the mid-infrared (A > 3/^m). Currently the most popu- 
lar model of the infrared (IR) spectrum in QSOs is ther- 
mal emission from hot dust. Evidence includes the lack 
of IR variability in most QSOs and the ubiquitous spec- 
tral dip around one micron which can be ascribed to dust 
sublimation at a temperature of Ts ~ 1800K. However, 
some radio quiet QSOs show phased IR/optical variability 
(Neugebauer & Matthews 1999), and dust absorption fea- 
tures are not generally seen in the mid-IR in bright AGN 
(Roche et al. 1991). Alternatively, the IR spectrum could 
be dominated by nonthermal synchrotron emission from 
the region near the black hole. 

Microlensing variability in the mid-IR can place a limit 
on the size of the source, distinguishing between dust and 
synchrotron models. The minimum radius at which dust 
can exist near the central engine is 
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where L is the QSO luminosity, e is the dust emissivity, 
and Ts is the dust sublimation temperature. i?dust is larger 
than i?E so that if dust is responsible for the emission, 
microlensing fluctuations will be very small (~ 0.05 mag- 
nitudes for -Rdust/-RE — 30, Refsdal and Stabell 1991). 

On the other hand, if the infrared emission is due to 
synchrotron, then a turnover should be present at low fre- 
quencies due to self-absorption. Such turnovers are typi- 
cally observed in quasars in the far infrared, but they are 
much steeper than the Fi, oc v^^"^ behavior expected from 
homogeneous self-absorbed synchrotron models (Hughes 
et al. 1993). Nevertheless, if such a turnover is due to 
self-absorption, then the requirement that the source be 
optically thick at the turnover frequency places a strong 
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upper bound on the size of the emission region of 

«--'°"(T^)"^(i)'"(To^)'""- 

(2) 

where is the flux at the turnover wavelength A, and 
B is the magnetic field in the emission region. Unfortu- 
nately the Einstein Cross lies in a region of the sky which 
was not surveyed by IRAS, so the far infrared flux and 
turnover wavelength are unknown. The typical flux in the 
far-infrared for a QSO as bright as the Einstein Cross is 
~ 10 mjy, so the size of the synchrotron emission region 
is expected to be comparable to the optical size (where B 
could be ~ 10^ G). The limit on the size of a synchrotron 
emitting region is therefore much smaller than i?E- If syn- 
chrotron emission is responsible for the IR, then microlens- 
ing variations in this region of the spectrum will be at least 
as large as in the optical. 

A lack of microlensing might allow measurement of the 
amplification due to the gravitational lens galaxy. If the 
IR emission is due to dust, then it will not be microlensed, 
and, as emphasized by Kochanck (1991), would provide 
a direct measurement of the relative magnifications due 
to the lens galaxy rather than due to microlensing (al- 
though microlensing by clusters of stars can still affect a 
large emission region). This data is required for accurate 
modeling of this lens (Witt and Mao 1994). This has been 
attempted by using the CIII] line, which originates from 
the extended broad line region and is therefore not mi- 
crolensed (Yee & De Robertis 1992, Racine 1992, Fitte 
& Adam 1994, Saust 1994, Lewis et al. 1998); however, 
the CIII] line is affected by extinction and continuum sub- 
traction. To avoid these problems, the flux ratios were 
measured in the radio by Falco et al. (1996) with a signal- 
to-noise ratio (SNR) of only 2-4. One can do much bet- 
ter in the IR. Assuming a power-law mass distribution in 
the lensing galaxy (M ex r^), Wambsganss & Paczyhski 
(1994) demonstrated that the total magniflcation of the 
lens could vary dramatically as f3 varies. Their model pre- 
dicts different flux ratios between the four images as f3 
varies, so measuring the flux ratios of an extended source 
might constrain the total ampliflcation. 

We present the first observations of the Einstein Cross 
near 10 microns (~3 microns rest frame). These obser- 
vations were taken at Keck, which is the only telescope 
with a large enough collecting area to detect such a faint 
QSO. In addition, at 10 microns it is diffraction limited 
with about 0'.'25 resolution, making it easy to distinguish 
between the four images which have a smallest separation 
of about |l[.' In section 2, the observations are discussed. In 
section 3, the data reduction is discussed and the results 
are presented. In section 4 the implications for emission 
models and lens models are discussed. In section 5 the 
results are summarized and future directions are pointed 
out. 

2. OBSERVATIONS 

The Einstein Cross was observed with the Keck Long 
Wavelength Spectrometer (LWS) on three nights in 1999. 
The LWS is a 128 x 128 Boeing Si: As array with a 10'.'5 
field of view at the Cassegrain focus of Keck I (Jones & 
Puetter 1993). The LWS was used in imaging mode with 
chop and nod which were necessary due to the faintness 



of the source and sky variability. The chop and nod were 
both fixed at 10", which is off of the LWS chip. The field 
was dithered by 1" in a 5-point pattern. The 8.9 and 
11.7 /im filters were used, which have ^ 80% flux trans- 
mission from 8.4 to 9.2 /im and 11.2 to 12.2 fim, respec- 
tively. Total observing times and on-source integration 
times are listed in Table H- In addition to observing the 
QSO, the stars HR 5616, 4 Lac, and HR 8551 were ob- 
served (near the same airmass) at 11.7 microns, yielding 
13,300 ADU/Jy/sec with a standard deviation of 5% for 
different objects and different nights. We observed stars 
HR 8551, a Boo, and HR 2443 at 8.9 microns, yielding 
10,300 ADU/Jy/sec with a standard deviation of 6%. 




Fig. 1: 5" X 5" image of the Einstein Cross at 8.9 
and 11.7 microns taken with LWS on July 28, 1999 and 
September 25 & 26 1999 (all data summed together; North 
is up). The images are labeled according to Yee (1988). 
bottom-right image. 



3. REDUCTION 

3.1. Basic Reduction 

Successive chopped frames were subtracted and added 
to successive nod sets to get each frame. Frames were me- 
dian clipped to get rid of cosmic rays and bad pixels. As 
the QSO is rather faint, it is not detected in any single 
frame, which typically consists of about flve minutes of in- 
tegration. To detect the QSO, the data were de-dithered 
and successive frames were added. Due to imperfect point- 
ing, the de-dithering can cause some broadening of the 
point-spread- function (PSF). As the QSO is too faint to al- 
low correction of pointing errors for individual frames, the 
stellar PSF was convolved with gaussians of various widths 
to see how this affects the results, as discussed below. The 
standard star PSF was constructed using an identical de- 
dithering procedure as that applied to the QSO images, 
but also was corrected for pointing errors. Scaled and 
broadened stellar PSFs were fltted and subtracted using 



Agol, Jones, & Blaes 



3 



Table 1 
Observing times 



Date 


Filter 


Observing time 


On-source time 


7/28/99 


11.7 fim 


50 min 


15 min 


7/28/99 


8.9 Aim 


43 min 


12 min 


9/25/99 


11.7 fim 


157 min 


45 min 


9/26/99 


8.9 Aim 


60 min 


16 min 



the known relative QSO image positions (from Falco et al. 
1996) while varying absolute position on the chip to min- 
imize the residuals of the image minus the model using a 
linear least-squares fit. Using the relative positions mea- 
sured by HST (Rix et al. 1992) changed the flux ratios by 
an amount much smaller than the error bars. The resulting 
fluxes are shown in Table ^j. Figure 1 shows the summed 
image from both bands and both nights, smoothed with a 
0.25" boxcar to eliminate noise. The lens galaxy is unde- 
tected; neither is there any evidence of a fifth lensed QSO 
image. 

For comparison to models, the fraction of the total Ein- 
stein Cross flux contained in each image was computed, 
hereafter referred to as the "flux ratio" ; i.e. if all four im- 
ages have the same flux, then the flux ratio is 1/4 for each 
image. 

3.2. PSF 

As mentioned above, the pointing errors resulting from 
dithering can lead to a broadening of the PSF for long ob- 
servations on faint sources, such as those presented here. 
A stellar PSF (corrected for pointing errors) was convolved 
with a gaussian, varying the FWHM until the best fit was 
obtained for each image. The typical FWHM varies be- 
tween 3 and 5 pixels (see Table ||) , but the flux ratio of each 
image on each night in each waveband varies by only 5% 
if the FWHM is changed from 3 to 5 pixels (0'.'25 to 0'.'4). 
This demonstrates that the relative flux ratios are rather 
insensitive to the gaussian FWHM. The absolute fluxes, 
however, are quite sensitive to the assumed FWHM. For- 
tunately, the argument for microlensing depends only on 
the relative flux ratios. 

3.3. Error bars 

A Monte Carlo technique was used to compute the er- 
ror bars on the fluxes. During these observations the LWS 
chip exhibited pattern noise which prevented a direct mea- 
surement of the error bars from photon counting statistics. 
Instead, the best-fit model of the four images was added 
back in at various points on the chip to create simulated 
images. The simulated images were run through the en- 
tire reduction procedure, measuring the position, gaussian 
FWHM, and fluxes of the four images. This was repeated 
~100 times to compute the standard deviation of the im- 
age fluxes, flux ratios, FWHM, and position. The error 
bars on the flux ratios are typically about a factor of 2 
greater than what one would infer from shot-noise statis- 
tics only. Thus, only the Monte Carlo error bars are re- 
ported in this paper (Table 0; 1-sigma error bars are al- 



ways quoted). The error bars on the absolute flux are quite 
large, mostly due to the uncertainty in the point spread 
function. The error bars on the measured FWHM are re- 
ported in Table ||; this is the FWHM gaussian which was 
convolved with the stellar PSF to approximate the QSO 
PSF. 

3.4. Consistency of the data 

The flux ratios measured on both nights and in both 
bands are consistent with being constant with time and 
with wavelength, as shown in Figure 2. Any disagreement 
between the measurements are consistent with random er- 
rors causing the differences; there is no need to invoke 
variability in the flux ratios with time or between the two 
bands, nor any need to find systematic errors to explain 
the differences. 
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Fig. 2: Fraction of the total flux of the Einstein Cross 
contained in each image. The data from September are 
plotted in bold; the 11.7 /im data with asterisks; and the 
8.9 /im data with squares. The error bars reflect the sta- 
tistical uncertainties 



4. DISCUSSION 



4.1. 



Microlensing Constraints 

Table I compares the mid-infrared flux ratios (averaged 
from both nights and both bands) to the flux ratios in 
other wave bands. The V band magnitudes are from the 
OGLE monitoring data on the dates closest to the Keck 
observations (available on their website; Wozniak et al. 



4 Mid-IR Imaging of QSO 2237+0305 



Table 2 
Infrared fluxes 





Flux 


Flux 


Flux 


Flux 


Image 


8.9/im 


11.7/im 


8.9/im 


11.7/im 




(mJy) 


(mJy) 


(mJy) 


(mJy) 


Date 


7/28/99 


7/28/99 


9/24/99 


9/24/99 


FWHM (pixels) 


4.5±0.6 


3.1±0.8 


4.5±0.6 


4.9±0.6 


A 


6.4±1.2 


5.7±1.0 


6.5±1.4 


6.4±1.5 


B 


9.1±1.2 


6.1±0.9 


5.7±1.3 


5.8±1.3 


C 


4.2±1.0 


3.5±0.7 


3.3±1.2 


2.6±1.2 


D 


7.0±1.1 


6.2±0.9 


6.8±1.4 


4.2±1.4 


All 


26.8±2.7 


20.5±2.1 


22.3±3.3 


18.9±3.6 



2000); this data was corrected for extinction as described 
in the Appendix. The VLA radio data are from Falco et 
al. (1996) which were taken in 1995. The GUI] ratios are 
from Racine (1992), data taken with CFHT in 1991. One 
can immediately draw the conclusion that the infrared is 
extended (and thus probably dust emission) from the fol- 
lowing lines of evidence: 

1) The flux ratios diff'er between the IR & V band by 
more than 6a for each image. Simultaneous with the Keck 
observations were strong optical microlensing events in im- 
ages A and C at the time of the infrared observations 
(OGLE website; Wozniak et al. 2000). Since the opti- 
cal is known to be microlensed, this indicates that the IR 
is extended on a scale comparable to or larger than an 
Einstein radius, causing dramatic differences in the flux 
ratios. 

2) The image ratios in the IR are close to the ratios 
predicted by some lens models, lending weight to the 
idea that the IR is more extended and thus not affected 
by microlensing. Table |^ shows the predictions of the 
A = 1 model from Schmidt, Webster, & Lewis (1998, here- 
after SWL), which were derived without using the flux 
ratios as a model constraint. For all published models 

= J2i={A,B,c,D}iMo^^)- M'^'^d,))'^ /'^i was computed, 
where fi is fractional flux of each image and ct; is the mea- 
sured error bar. The three best models are A = 1 of SWL 
with = 7, model 4 of Kochanek (1991) with — 7, and 
model 1 of Witt, Mao, & Schechter (1995) with = 8. 
These x^ indicate that there may be deviations of the in- 
frared magnifications from the predicted magnifications at 
the level of 10%. This may be due to inaccuracies in the 
modeling, an underestimation of the error bars, or mi- 
crolensing fluctuations at the 10% level, which is to be 
expected for a source 15 times the size of an Einstein ra- 
dius (Refsdal & Stabell 1991, assuming a — 0.5). All other 
published gravitational lens models for Einstein Cross have 
larger (Kent & Falco 1988; Kochanek 1991; Rix, Schnei- 
der, & Bahcall 1992; Wambsganss & Paczyhski 1994; Witt 
et al. 1995; Chae, Turnshek, & Khersonsky 1998). 

3) The image ratios at 8.9 and 11.7 microns are consis- 
tent within the error bars, indicating that neither wave- 
length is affected by either extinction or differing mi- 
crolensing amplification due to different sizes as a function 



of wavelength. 

4) The radio emission is thought to be extended on a 
scale larger than an Einstein radius (possibly in a jet) and 
thus is unaffected by microlensing. Since the IR image 
ratios are consistent with the radio image ratios within 
the error bars, this implies that the IR may also be un- 
affected by microlensing. Note that the infrared SNR is 
much better than in the radio. 

5) In the case of images C and D, the IR image ratios 
are not consistent with the GUI] image ratios, but agree 
better than with the V band flux ratios. The deviations 
may be due to microlensing fluctuations of GUI] if it is 
only a few times larger than an Einstein radius, or due to 
the systematic errors which affect the GUI] measurements. 

6) Image G dropped in optical flux by a factor of 1.4 
between August 4 and September 26 (the two OGLE ob- 
servations closest in time to the two Keck observations). 
Our observations are consistent with no variation in im- 
age G in either band between the two dates, although the 
possibility of a variation is not strongly ruled out. 

7) Gonsider a model in which a point source and an ex- 
tended source contribute to the infrared flux (which would 
be the case if both a dusty torus and synchrotron source 
contributed to the infrared flux). Then, the synchrotron 
source would be magnifled by the same amount as the 
optical (due to microlensing), while the extended source 
would be magnifled by the macrolensing magniflcation (as- 
suming the SWL model holds). With these assumptions 
one can compute what the fraction of point source flux is 
relative to the extended source flux in order to give the 
infrared flux ratios. Carrying out this exercise, one flnds 
that the point source contribution must be negative to ex- 
plain the difference between images A and B or A and G. 
Simply put, the infrared flux ratio of B is larger than the 
model, while the point source (i.e. optical) is smaller than 
the model (see Table ||), so adding a point source to an 
extended source can only decrease the flux with respect 
to the model. Since the point source flux must be non- 
negative, it is most likely zero. Then, the disagreement 
between the infrared flux ratios and the model flux ratios 
is probably due to inaccuracies in the lens model, Poisson 
fluctuations due to a random number of lenses passing in 
front of a given image (since the source is not truly infl- 
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nite), and/or errors in the measured fluxes. 

Given these lines of evidenee, one conehides that the 
mid-infrared emission region is larger than the Einstein ra- 
dius of a typical star in the lens galaxy, about 10^^ cm. The 
infrared may be extended by 0'.'24 which corresponds to 
about 0.8 kpc in the source plane; however, this extension 
is possibly due to pointing errors as discussed above. Con- 
sequently, 0.03pc ^ RiR <^ 800pc. Comparing the optical 
and infrared data with simulated microlensing lightcurves 
with various source sizes, one can compute the probability 
distribution of the infrared source size, which will be done 
in future work (Wyithe & Agol, in preparation). 




log(KHr)) 



Fig. 3: Spectral energy distribution of the Einstein 
Cross (non-simultaneous; frequency in observed frame). 
The near-infrared and optical points are taken from 

Nadeau et al. (1991); the ultraviolet from Blanton, 

Turner, & Wambsganss (1998); the X-ray data (dashed 

line) is from Wambsganss et al. (1999); the radio points 

are from Falco et al. (1996); the mm points (diamonds) are 

3 CT upper limits from Richard Barvainis (priv. comm.); 

and the mid-infrared data are from this work. The solid 
line is the radio-quiet QSO composite from Elvis et al. 

(1994) scaled to match the optical flux. The dotted line is 
the dust shell model described in the text based on Bar- 
vainis (1987). 

4.2. Spectral Energy Distribution 

To study in more detail the properties of the Einstein 
Cross, an SED was constructed of all available data from 
the radio to X-ray from VLA, IRAM, CFHT, Keck, EST, 
and ROSAT. The lens galaxy has been removed (or does 
not contribute) for all the observations shown. The near- 
IR, optical, UV, and X-ray were corrected for extinction 
in the lens galaxy and the Milky Way as described in the 
Appendix. This brings the optical/UV slope into rough 



agreement with the average QSO slope. Figure 3 shows 
the SED in the observed frame with the extinction correc- 
tions, but without a correction for magnification. Figure 
3 also shows the Elvis et al. (1994) radio-quiet QSO com- 
posite redshifted and scaled to the R band flux. Assuming 
a magnification of 16 (from the SWL model - this value is 
used throughout the rest of this paper), this implies a QSO 
luminosity of 2 x 10*^ erg/s for the scaled composite. How- 
ever, there appear to be deviations of the Einstein Cross 
from the composite QSO in the mm, infrared, UV, and 
X-ray. The mm upper limits are slightly lower than the 
average SED. The infrared has a narrower peak at higher 
frequency than the average QSO. A similar peak exists in 
3C 273; however, 3C 273 peaks at 3 /zm and is broader 
than the Einstein Cross (Barvainis 1987). The ultraviolet 
points are high, which may mean that this is a hot QSO or 
that the QSO has varied (or been magnified by microlens- 
ing) as the UV points were taken nearly a decade after the 
optical/near-IR points. The X-rays fall a factor of 3 be- 
low the average QSO, which is within the large scatter of 
optical to X-ray flux ratios of QSOs. The radio to optical 
flux ratio is identical to the QSO composite. The infrared 
discrepancy is of greatest relevance to this paper, which is 
discussed in turn. 

4.3. Origin of Infrared Bump 

The Keck mid-IR fluxes (summed over all four im- 
ages) give a spectral index of a = 0.5 ± 0.2 [F^, oc 
which is quite different from the average for PC QSOs, 
a = -1.4 ±0.5 (Neugebauer et al. 1987). This unusual 
behavior may be due to calibration problems. However, 
we think our error bar estimates are accurate, and we are 
unable to identify any calibration errors which affect the 
Keck observations at the 50% level which would be re- 
quired to bring the slope into agreement with that of an 
average QSO. The data were compared with ISO CAM ob- 
servations (TDT 54500601; P.I. Sangeeta Malhotra) of the 
Einstein Cross from the ISO archive which show a similar 
infrared slope of a = 0.45. The ISO CAM fluxes at 9.1 and 
11.7 microns agree with the LWS measured fluxes within 
the 1(7 error bars. The near-IR spectral index is dramat- 
ically different: Nadeau et al. (1991) report a = —3.7 
between the K-band and L-band (3.3 micron band). This 
implies that there is a peak in i/F^ between 3.3 microns 
and 8.9 microns (1.2 microns and 3.3 microns in the rest 
frame); i.e. the peak infrared flux occurs somewhere near 
2 microns. Typically, QSOs peak much further in the in- 
frared, near 10 microns (see Figure 3). 

This unusual infrared bump might be explained as a 
quasi-blackbody component. There are several candidate 
physical components which might contribute near these 
wavelengths: synchrotron, hot dust, a star burst, or the 
host galaxy. Synchrotron is ruled out by microlensing, so 
the other possibilities are considered next. 

4.3.1. Hot Dust 

A blackbody component was flt at the dust sublimation 
temperature (~ 2000 K), which agrees with the data re- 
markably well if the effective radius (area = 7rr^) of the 
source is = 2 x 10^® cm. Treating the dust emission as 
a blackbody assumes the emissivity is unity, appropriate 
for large grains; small grains would reduce the emissiv- 
ity, requiring a larger area, but lowering the temperature 
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Table 3 
Flux fractions of images 



Wavelength 


Date 


Image A 


Image B 


Image G 


Image D 


8.9 & 11.7^m 


7-28 & 9-24 


0.27±0.02 


0.30±0.02 


0.16±0.02 


0.27±0.02 


V-band'' 


8-4-99 


0.39 


0.11 


0.41 


0.10 


V-band^ 


9-26-99 


0.46 


0.12 


0.32 


0.10 


8 GHz 


1995 


0.29±0.06 


0.32±0.16 


0.16±0.06 


0.22±0.06 


GUI] 


1991 


0.26±0.01 


0.29±0.01 


0.26±0.01 


0.19±0.01 


Model (SWL) 




0.25 


0.27 


0.15 


0.32 



"^Corrected for extinction 



to 1600 K (for e = 27rr/A, where r ^ Ijim is the dust 
grain radius). The inferred source size agrees with the 
dust sublimation estimate and with the microlensing size 
constraint. However, this implies that the dust should be 
located at a constant radius so that when it is exposed to 
the QSO continuum it equilibrates at a constant temper- 
ature. This means that the dust is in a shell geometry; if 
there is dust at a range of radii, then it will radiate at lower 
temperatures than the sublimation temperature which will 
make the infrared SED wider and flatter than a blackbody, 
typically rising in flux towards longer wavelengths (Pier & 
Krolik i992). 

To quantify the constraint on geometry, the near- 
infrared bump was fit with the spectral model of Barvainis 

(1987). Barvainis' model assumes that the dust is in a shell 
subtending fl steradians of the central source with an in- 
ner radius (ri) set by the luminosity of the source (L) and 
the dust sublimation temperature (Tmax)- The shell is as- 
sumed to be optically thin in the ncar-IR (A ^ 1m™)j with 
a UV optical depth of 3. Barvainis' imiform dust model 
was adopted with a = 0.05/xm, emissivity proportional to 
u^'"^, and we additionally assume that the dust grains have 
a constant number density, that L = 2 x IC*^ erg/s, and 
that Tmax — 1800K, fixing the inner radius at l.l pc, con- 
sistent with the microlensing size constraint. The two free 
parameters left are r2/ri (the ratio of the outer to inner 
radius of the shell) and fl. An adequate fit was found for 
r2/ri = 3 and fl = Ttt/A. In this model, the dust outside 
r2 is shielded, and so does not contribute to the repro- 
cessed light. The grain number density at the inner edge 
is 5 X 10^ cm~^, so the implied dust mass is about 4M0 
for a dust density of 1 g/cm^. The Toomre Q-parameter 
is quite small, ^ 10^^ (for a mass of 2 x 10^ Mq), which 
indicates that this region should be gravitationally unsta- 
ble. A clumpy disk model should have a similar spectrum 
since the smaller filling factor can be compensated by an 
increase in dust number density. The conclusion is that the 
narrow infrared peak requires a narrow range of radii (fac- 
tor of 3) so that reprocessed luminosity is absorbed by a 
region with a narrow range of dust temperature. This nar- 
row range is not unprecedented as Barvainis (1992) found 
a narrow range of radii for the dust emission in Fairall 9 
from reverberation ri <^ 0.3pc and r2 ~ 1.3pc. 

This, of course, is a rather simplistic model. X-rays may 
penetrate deeper, heating dust at larger radii. Transient 



heating by X-rays causes temperature flickering which 
broadens the spectrum; this is only important for dust 
grains with radii a ^ 0.01/Ltm (Voit 1991). Infrared radi- 
ation transfer has been ignored since the shell is optically 
thin by construction; reality may not be so simple. A 
warped disk with a warp angle which scales logarithmi- 
cally radius would produce a much wider infrared bump 
than observed. 

4.3.2. Host Galaxy or Starhurst 

Normal galaxies can be well described by a 3500 K black- 
body in the near infrared spectral region (Schmitt et al. 
1997), which peaks at a wavelength of about Ijim (rest), 
while the Einstein Gross spectrum peaks closer to 2/zm 
(rest). Thus, to fit the infrared bump requires strong ex- 
tinction as in a starburst galaxy: an extinction of Ay ~ 7 
is needed to fit the infrared bump; most of the light is ab- 
sorbed by dust. The absorbed luminosity is Labs = 6 x 10**^ 
erg/s, which must be re-radiated in the far-infrared (FIR). 
The radio-FIR relation (Gondon 1992) limits the FIR lu- 
minosity to 6 X lO"'^ erg/s, a factor of 10 times smaller 
than the absorbed luminosity. A test of the host galaxy 
hypothesis would be to look for an extended source in the 
IR images. Meurer et al. (1997) have found an empiri- 
cal limit on the bolomctric surfacx; brightness of starburst 
galaxies of 2 x lO^^L0/kpc^, which would indicate that this 
starburst has a size of ^ 3.5 kpc, or ~ 1" (including the 
factor of ^ 2 stretch due to Icnsing). This is larger than 
the infrared upper Imit of 0'.'25, ruling out a starburst as 
the origin of the infrared bump. 

One can roughly estimate the host galaxy luminosity 
from the QSO luminosity as follows. Ranch & Bland- 
ford (1991) have argued that sub-Eddington accretion disk 
models are too large to be consistent with the optical size 
constraint from microlensing. We confirmed this by fitting 
the optical/UV SED with various thin accretion disk mod- 
els from Hubeny et al. (2000); all disks are a factor of a few 
too big to be consistent with the microlensing constraint. 
A smaller emission region might be achieved by decreasing 
the mass of the black hole; however, once the accretion rate 
reaches Eddington, the luminosity will be roughly fixed at 
the Eddington rate. Thus, if one assumes that the ob- 
served (isotropic) luminosity equals the Eddington rate, 
then the black hole mass is about 2 x IO^Mq. Using the 
Magorrian et al. (1998) relation between black hole mass 
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and bulge luminosity, the bulge luminosity is estimated to 
be about 6 x IO^Lq, with an uncertainty of about a factor 
of 10. This is more than 150 times smaller than the near-IR 
luminosity, yet another indication that the infrared bump 
may not be due to a quiescent host galaxy. 

5. CONCLUSIONS 

The observations presented here indicate the mid- 
infrared source size in the the Einstein Cross must be 
greater than lO^^cm, ruling out synchrotron emission. If 
true for other radio-quiet QSOs, this implies that fits to 
the Big Blue Bump which include a synchrotron power law 
may be flawed (Sun & Malkan 1989, Laor 1990). Starburst 
and host galaxy contributions to the mid-infrared emission 
are ruled out based on the low radio flux, small size, and 
large infrared luminosity, leaving hot dust as the only vi- 
able model for the infrared emission in the Einstein Cross. 
To avoid sublimation, the dust must be at a distance of 
about 1 pc from the QSO; to explain the narrowness of 
the infrared bump it must be located within ~ 3pc (ac- 
cording to a simple reprocessing model); and to explain its 
strength, the dust must intercept about half of the QSO 
continuum. Since the mid-IR is extended, the flux ratios 
should be nearly equal to those of the galaxy gravitational 
lens model, so future models should be constructed with 
the infrared flux ratios as model constraints. Improved mi- 
crolensing statistics, magnifications, and time delays will 
result from better lens models. Accurate monitoring of 
the Einstein Cross at mid-infrared wavelengths with, say, 



NGST, may reveal flux variations at the few % level, which 
would allow a better estimate of the size of the dust emis- 
sion region. The Einstein Cross mid-infrared spectrum is 
quite different from any other QSO; it should be observed 
with SIRTF to better constrain the dust emission model. 
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APPENDIX 

APPENDIX; CORRECTION FOR EXTINCTION 

To correct for extinction in the lens galaxy, the correlation between g — i color (Gunn bands) and lens galaxy surface 
brightness found by Racine (1991) was used. Extrapolation to zero surface brightness then gives the color corrected for 
extinction in the lens galaxy, if one assumes that the dust column is proportional to surface brightness (this assumption 
may be invalid if the dust is patchy). To solve for the absolute extinction, we assume the extinction curve is described 
by i? = 3.1 ± 0.9, the average Milky Way extinction curve. Ay for image A is 0.88±0.21, B is 0.84±0.20, C is 1.30±0.31, 
and D is 1.15±0.27. Using the Milky Way extinction E{B — V) = 0.071 measured by Schlegel, Finkbeiner, & Davis 
(1998), the Nadeau et al. spectrum was corrected for each image for galactic extinction using the R = 3.1 extinction 
curve (Fitzpatrick 1998). It is reassuring that the corrected spectrum is similar to that of an average QSO (see Figure 3). 
Yee (1988) and Falco et al. (1999) correct for extinction assuming the Einstein Cross has a spectral index equal to the 
average value for QSOs; since individual QSOs do not behave as the average, we prefer the Racine technique. 

To correct the X-ray data for absorption, we adopt the dust-to-gas ratio given by Spitzer (1978): Nh = 5.9 x 
10^'^£'s-\/mag~^cm~^. An average extinction of 1 for the lens galaxy and 0.2 for the Milky Way implies Nh = 
2.5 X lO^^cm^^. The X-ray power law was assumed to be equal to that of an average QSO and the count rate of 
Wambsganss et al. (1999) was used to compute the unabsorbed X-ray flux in the ROSAT band. 
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